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The production of particles with double, triple and hidden charm in heavy ion collisions is studied 
in the framework of the statistical coalescence model. According to the postulates of the model, the 
charm quark-antiquark pairs are created at the initial stage of a heavy ion reaction in hard parton 
collisions. The amount of charm is assumed to be unchanged at later stages. The charm (anti)quarks 
are distributed among different hadron species at hadronization according to the laws of statistical 
physics. Several approaches to the statistical treatment of charm hadronization are considered. The 
grand canonical approach is appropriate for systems containing large number of charm (anti)quarks. 
The exact charm conservation and Poissonian fluctuations of the number of charm quark-antiquark 
pairs should be taken into account, if the average number of these pairs is of oder of unity or smaller. 
The charm hadronization in a subsystem of a larger system is discussed. It is explained why the 
canonical approach is not appropriate for the description of charm hadronization. The obtained 
formulas can be used to calculate the production of charm in heavy ion collisions in a wide energy 
range. 
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I. INTRODUCTION 

The thermal hadron gas (HG) model has demonstrated 
an obvious success in describing the chemical composi- 
tion of light-flavored hadrons produced in heavy ion |l| 
and even in elementary hadron-hadron as well as in 
electron-positron Q collisions. The experimental data 
can be well fitted with only three free parameters: the 
temperature T, the volume V and the baryonic chemical 
potential ^ of the hadron gas at the point of the chemical 
freeze-out. (Sometimes the fit is improved by introducing 
one more parameter — the strangeness suppression fac- 
tor 7 S .) This success motivated attempts to extend the 
applicability domain of the thermal model also to heavy- 
flavored hadrons, for instance, to the description of J/ip 
meson production 

A straightforward application of the equilibrium HG 
model to hadrons with open and hidden charm is not, 
however, justified. Partons with rather large momenta 
are needed to produce a heavy quark-antiquark pair. 
Consequently, the time of charm equilibration in a ther- 
mal hadronic or even quark-gluon medium is large and 
definitely exceeds the lifetime of the fireball. Produc- 
tion of charm can take place only at the initial stage 
of the heavy ion reaction, when hard partons are avail- 
able. The charm production/annihilation rate is too low 
to keep the number of heavy quark-antiquark pairs at 
its chemical equilibrium value at later stages. There- 
fore, the total amount of heavy flavor should be out of 
equilibrium at the point of hadronization and chemical 
freeze-out. It is reasonable to expect, however, the dis- 
tribution of heavy quarks and antiquarks among different 
hadrons with open and hidden charm to be thermal and 
controlled by the same values of thermodynamic param- 
eters that fit the chemical composition of light-flavored 



hadrons. These ideas are implemented in the statistical 
coalescence model (SCM) UllQ- 

In the present paper, I consider SCM, which is based 
on the following postulates: 

• The charm quarks, c, and antiquarks, c, are created 
at the initial stage of A+A reaction in hard parton 
collisions. 

• Creation and annihilation of cc pairs at later stages 
can be neglected. 

• The formation of observed hadrons with open and 
hidden charm takes place near the point of chemical 
freeze-out in accordance with the laws of statistical 
physics. 

This approach appeared to be quite successful describ- 
ing the J/ip and ip' production in (semi)central Pb+Pb 
collisions at Super Proton Synchrotron (SPS) Al- 
though the role of the statistical coalescence at SPS ener- 
gies is still under discussion its dominance at the en- 
ergies of the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC) is more certain. The 
SCM-based predictions for the quarkonium production at 
RHIC [lllllj are consistent with transport calculations 


Quarkonia are not the only type of hadrons whose pro- 
duction can be described by the statistical coalescence 
model. If two, three or more charm quark-antiquark pairs 
arc created in a nucleus- nucleus collision, hadrons with 
double and triple charm can be formed. These hadronic 
states were long time ago predicted by the quark model 
(see and references therein), but most of them have 
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not been observed experimentally yet. 1 

The hadrons with double and triple charm are of 
special interest from both theoretical end experimental 
point of view. Due to a rather large mass of charm 
(anti)quarks, their interactions within a hadron are close 
to the perturbative regime of Quantum Cromodynamics 
(QCD). A study of the properties of doubly and espe- 
cially triply charmed baryons would allow to test QCD- 
based models of quark forces. An observation of more 
exotic hadronic states, like multi-charm tetra- and pen- 
taquarks will open a new window into the structure of 
the hadronic matter. The intention to detect the double 
and triple charm poses a new challenge to the experi- 
mentalists and will demand a further development of the 
experimental technique. 

In the present paper, I derive the formulas that allows 
to calculate the yield of double and triple charm particles. 
The formulas for the hidden charm are also included as 
their detailed derivation has not been published yet. 

The article is organized as follows. In the Section [n] 
I consider the grand canonical approach. A system with 
exactly fixed numbers of charm quarks and antiquarks is 
studied in Section lTTTI A more realistic situation, the sys- 
tem with Poissonian fluctuations of the number of charm 
quark-antiquark pairs, is considered in Section llVI Sec- 
tion is devoted to the case, when only a part of the 
total system is available for the observation: the double, 
triple and hidden charm yield in a subsystem is studied. 
Summary is given in Section IVllI 



II. THE GRAND CANONICAL APPROACH TO 
THE STATISTICAL COALESCENCE 

The grand canonical version of the statistical coales- 
cence model was proposed in [5| ■ It can be applied to the 
systems containing large (A C c 3> 1) number of cc pairs. 

Let us start from the grand canonical partition func- 
tion for the ideal hadron gas in the Boltzmann approxi- 
mation: 

z(v, t, { A}) = e E E n [ lH 7; m) ] . 

where L is the total number of hadron species (including 
resonances), V and T are the volume and temperature, 
4>(T] m, g) is the one-particle partition function 

J>(T;m,g) = ^ j™ p*dp exp (- V ^ + ^ (2) 



1 The recent observatio n of the doubly charmed baryon SJ^, ll ■"! is 
still under discussion llfl . 



Here m is the particle mass, g is the degeneration factor 
(the number of spin states), K 2 is the modified Bessel 
function. In the nonrelativistic limit m >> T the ex- 
pression takes the form 

<f>(T;m,g) ~ sf-^r] exp(-m/T). (3) 

The fugacity A/ is expressed via the chemical potentials 
(electric — )i e , baryonic — /ib-, strangeness — fi s and 
charm — /i c ) and suppression (enhancement) factors (for 
strangeness — 7 S , and charm — 7 C ): 

A, = T M< 7 N< ex p(^), (4) 
[ii = qin e + bi^b + sifM s + cni c , (5) 

where qi, bi, sz, q, are, respectively, the electric charge, 
baryon number, strangeness and charm of the hadron 
species I, \s\i and \c\i are the numbers of valence strange 
and charmed (anti-)quarks. The chemical potentials in 
the right-hand-side of (0 are responsible for keeping the 
correct average values of the corresponding charges in 
the system. The suppression (enhancement) factors 7 S 
and 7 C are introduced to take into account a deviation 
of the total number of strangeness and charm from their 
equilibrium values 2 . 

The average number of particles is given for each 
species as 

dlogZ(V,T,{\}) 
Ni=M — = Xi(j){T,mi,gi)V. (6) 

The total number of hadrons 

Ntot — Ni, (I runs over all hadron species), (7) 

can be broken up into several pieces: ./Vo, the number of 
zero charm hadrons (excluding hidden charm), Nh, the 
number of hidden charm mesons, N±, Aj, N 2 , Aj, A 3 , 
Ag, the numbers of hadrons with, respectively, single, 
double and triple charm and anticharm: 

Ntot = N + N H + Ai + Ni + A 2 + A2 + A 3 + A 3 . (8) 

Let us consider a system containing in average N c 
charm quark and A g charm antiquarks. Then, as far as 
charm creation and annihilation are neglected, the fol- 
lowing equalities should be satisfied: 

(A c ) = m + NH + 2A 2 + 3A 3 , (9) 
(A 5 ) = Ai + A ff + 2A2 + 3A g . (10) 



2 The suppression (enhancement) factor 7 S (the same for 7 C ) is 
equivalent to an additional chemical potential which, in con- 
trast to fi s , has the same influence on the strangeness as on the 
antistrangeness: 7 S = exp (—5^) 
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It is easy to see from 0,(0) and © that the above equa- 
tions can be rewritten as 

(JV C ) = 7c A c iV 1 + 7c 2 iVH + 27 c 2 A2iV 2 +37 c 3 A^ 3 , (11) 
(N s ) = icK 1 ^ + i 2 c Nh + 2 1 2 X^ 2 N 2 + 3 7c 3 A c ^ 3 7V3, 

(12) 



where 



- frV»)*N H =(N c )(N £ )-£^-, (21) 

( 7 (°)) 2 (A(°)) 2 7V 2 = (iV c ) 2 §, (22) 
Nf 

(7 (0) )2(A (0))2 #3 = {Nc) S^ (23) 



N H 



Nk = AT fe j Mc=0) , 



k = #,1,1,2,2,3,3, (13) 



and 



Xr 



exp 



V T 



(14) 



The constituent charm quark mass is about an order of 
magnitude larger than the typical temperature of chem- 
ical freeze-out. From this reason, 



JVi, iV T > N H , N 2 , N 2 > N 3 , N 3 



(15) 



due to the exponential factor in J2J). Therefore, if the 
factor 7c is not extraordinary large, the terms in 111|) 
and (|12|) corresponding to hidden, double and triple 
(anti)charm can be neglected in a zero approximation 
and the coupled equations (|11|) and (|12fl can be simpli- 
fied: 



<JV„> = 7 f A^iVr, 

m = 7 f (aW)- 1 ^. 

The solution can be easily found: 
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(N C )(N 5 ) 
NiNi 

(Nc)^ ' 



(16) 
(17) 



(18) 



(19) 



Now one can calculate the number of single, hidden, dou- 
ble and triple charm particles in the zero approximation: 



("A = #i = (N c ), 



(20) 



The formulas for anticharm can be obtained by the ob- 
vious replacement c — ► c, 1 — > 1 etc. 



The formulas (|2UI - I23|) are obtained under the assump- 
tion that the charm enhancement factor 7s is not large, 
so that only a tiny fraction of the charm quarks and an- 
tiquarks hadronizes into hidden, double and triple charm 
particles. This may be not true at very high energies. 
Indeed, let us consider for example the leading order ex- 
pression H22[) for the number of double charm particles. 
The fraction of charm quarks that hadronize into double 
charm particles is proportional to the number of charm 
quarks in the system and inversely proportional to the 
system volume at freeze-out: 



(N c ) V 



(24) 



The volume is proportional to the multiplicity of light 
hadrons. The charm production cross section grows with 
the collision energy faster than the multiplicity of light 
hadrons. Therefore, the ratio l|24|> becomes comparable 
to unity at some point. In this case the above approxi- 
mation does not work. 



The accuracy of the leading-order approximation l|20|) - 
(|23|l can be estimated by calculating the next-to-leading 

order corrections. Let us substitute 7c — » 7 c°' ) + 7 1^ and 



A c -> A c 0) + A c 1} into (JTTJ and O. Neglecting higher 



order terms like 7 e Ac. 1 ' 



l^N k (k 



#,2,2,3,3), etc. 
and taking into account (|T^)l . (|T7Jl and (|2*T]) (|2*5|l one gets 







h 3(iv 3 ; 


,(0) 
>GCE 




Nt 






(Nh^ce- 


H2(2V 5 >^- 


H 3(iv 3 : 


M 

>GCE 



(25) 

A (0)^^ A (0) 7 (l) _ 7 (0) A (l)^ = {Mh)gcE + 'W^GCE + 3{M3)GCE _ ( 2 g) 



These coupled linear equations can be easily solved: 



n) = _1 (0) ( (Nh) ( § ) ce + 2(N 2 )^ ce + 3(N 3 )^ ce (N H )^ CE + 2{N^° CE + 3(N- 3 )^ CE \ 
7c 2 7c I <iV c ) + (iV 5 ) i ' 1 j 



Now the next-to-leading order corrections to the number 
of particles with different charm content can be calcu- 
lated: 



results of the SCM will be studied in the subsequent sec- 
tions. 



({Nh)<§ } ce + 2{N 2 )<§l E + 3(N 3 )M E 

(29) 

( N i)gCE ) f3Ql 

(N e ) ) ' m 



(0) 2(N 1 )^ CE 
{N%)gce (N c ) ' 

(N) (o) HN^ce 
Wgce {Nc} ■ 



(31) 
(32) 



The higher order corrections can be obtained in a similar 
way. Still, the expressions become rather unwieldy. It is 
more reasonable to solve the coupled equations and 
(|12H numerically, if the next-to-leading order corrections 
(|2H1>— (E21) become too large. 

The above formulas allow to calculate the total num- 
ber of particles with given (anti)charm content (hidden, 
single, double or triple) . If the particle number of a single 
species is needed, it can be found in the following way. 
One can deduce from equations © and © that the 
number of particles of a single species I can be found from 



(33) 



where k = H, 1,1, 2,2, 3 or 3, depending on the charm 
content of the species I. 

The grand canonical approach is the simplest version of 
the statistical coalescence model. However, it does not 
reflect the real experimental situation sufficiently well. 
Indeed, we operated only with average values (N c ) and 
(N e ) in the above consideration. The correlations be- 
tween N c and N 5 were ignored. In reality, however, the 
quarks and antiquaries are produced in pairs. So that not 
only average numbers coincide (N c ) = (Nc), but also the 
numbers of quarks and antiquarks in every single event 
are equal: N c — N s . The influence of this fact on the 



III. SYSTEM WITH FIXED NUMBERS OF 
CHARM QUARKS AND ANTIQUARKS 

In this section, a system with fixed numbers of charm 
quarks and antiquarks is considered. This does not cor- 
respond exactly to what is observed in the experiment. 
Still a study of such a system would allow us to obtain 
an important intermediate result. Later, it will be used 
in more realistic calculations. 

Using Eq.©, one can rewrite the partition function 
(|T|) in a more compact form: 



oc oc oc L 



N! 1 



wt,{x}) = EE- £119- 



z 1= 0i 2 =0 i L = ll=0 



= ex P r£ N n = ex p ( N **) ■ ( 34 ) 

The total number of particles in the last expression can 
be represented according to JSJ, then 

Z(V,T,{\}) = exp(]TAU =JJ eX p(iV fc ), (35) 

V k I k 

k = 0,#, 1,1,2,2,3,3. 

Expanding all the exponents in the product, except the 
first one, into the Taylor series one gets 



Z(V,T,{\}) 



N! 



EE £11 9b (36) 



i^— (Hi— 23— 1 k 

k = #,1,1,2,2,3,3. 



Due to the property of The Kronecker delta 
^ m <5(m,rt) = 1, nothing changes if we multiply an ex- 
pression by a Kronecker delta and sum over one of its 
indices. Doing this twice on the expression (|36[) . one gets 



Nl k 



Z(V,T,{\}) = e N °J2 E ■ • ■ E E E lfi + 81 + 12 + *h + 'I + '2 + *3) II "I' ( 37 ) 



j H =0u=Q ig=l N c Ns 

k = #,1,1,2,2,3,3. 
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Then, after changing the summation order and using J2J, © and (|13|) . the above expression can be rewritten as 



(N c -N £ )^\ Z NcNs (V,T,{\}). 



(38) 



Here Zn c n £ ^ s the partition functions for the systems containing exactly N c charm quarks and antiquarks: 



OO OO OO 



N? 



in— Oil— io — 1 



(39) 



fc = if, 1,1, 2, 2, 3, 3. 



This function is canonical with respect to the exact conservation of the number of charm quarks and antiquarks and 
grand canonical with respect to the conservation in average of all other charges, whose values are controlled by the 
activities 



- Afc| M< _ 0i7c=1 . 

The number of summations in (|39|l can be reduced due to the Kronecker deltas: 



(40) 



rXis n h N 2 N i ^ 3 N i N 1 m 1 

z NcNs {v,T,{\}) = Y- M rY- 2 rY- 2 r Y- 1 rY . 3 , .\ .\ 



(41) 



iH=0 12=0 



i 2 =0 

r 



i 3 =0 " i 3 =0 



Here 



■max 
l H 


= min(N c ,Nc), 


(42) 




= [(N c -i H )/2], 


(43) 


if ax {i H ) 


= [(JV e -* H )/2], 


(44) 




= [(JV c -i H -2i a )/3], 


(45) 




= [(JV e -i H -2t a )/3], 


(46) 




= [N c - i H - 2i 2 - 3i 3 ], 


(47) 




= [Nc-i H -2k-Zk]- 


(48) 



(The square brackets mean here the integer part, i.e. [x] 
is the largest integer number that does not exceed x.) 

The average number of hidden (single, double, triple) 
(anti)charm particles in the system with fixed numbers 
of charm quarks and antiquarks is found as 



, CM; *-r 



~ dN k d\ogZ NcNs 



dXi 



dN k 



k = H, 1,1, 2, 2, 3, 3; (49) 
I runs over all hadron species of the type k. 



From 



one sees that 



therefore 



V ox. 



(50) 



(51) 



and finally 



(N h )fix = N k - 



1 dZ 



N c Nc 



J N c Nc 



ON, 



(52) 



I restrict my further consideration to the case, when 
the partition function l|41[l is dominated 3 by the term 
with in = ii = i 2 — *3 — *g — 0, i\ = N c and ij = N s : 



7 (o) 

J N C N E 



(53) 



iV c ! m ' 

The same term dominates also the derivatives of Zn c n e 
with respect to N c and N 5 . It is easy to see that in this 
case, most of the charm hadronizes into hadrons contain- 
ing only one c-quark or antiquark: 



(o) 



N c , 

Nr, 



(54) 
(55) 



Only a tiny fraction of the total charm is accommodated 
into hidden, double and triple charm hadrons. 

The leading term (|53|l does not depend on in- There- 
fore, the derivative dZM^N^/dNn is dominated by the 



term of (|41|l with in = 1, i<i 
i\ = N c — 1 and ij = N e — 1: 



13 



dZ 



(o) 

N c Nc 



N c -1 



N E -1 



ON, 



H 



(N c -iy.(N- c -iy. 



o. 



(56) 



3 Due to the inequality 1151 . the contribution of other terms be- 
comes sizable only at very large N c and/or Nc- This case can 
be treated within the grand canonical approach and there is no 
reason to study it here. 
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From (|52|) one finds the average number of hidden charm 
particles in the system with fixed numbers of charm 
quarks and antiquarks: 



NiNi 



(57) 



Similarly, the average numbers of double and triple 
charm particles are given by: 



N C (N C 



1)^4 



and 



iV c (7V c -l)(7V c 



N-i 
2)-^-, 



(58) 



(59) 



respectively. The corresponding formulas for anticharm 
can be obtained from (|58|l and l|59|) by the obvious re- 
placement c — > c, 1 — * 1 etc. 

We have considered the charm hadron production by 
quark coalescence in the thermal hadron system with 
fixed numbers of charm quarks and antiquarks. In real- 
ity, this number is not fixed. It fluctuates from event to 
event. These fluctuations have to be taken into account 
in more realistic calculations. 



IV. THE SYSTEM WITH POISSONIAN 
FLUCTUATIONS OF THE NUMBER OF CHARM 
QUARK-ANTIQUARK PAIRS 

In a relativistic collision of two nuclei, the charm 
quarks and antiquarks are created in pairs in indepen- 
dent nucleon-nucleon collisions. Therefore, the number 
of quarks in the system is always equal to the number of 
antiquarks: 



N r = Nr 



(60) 



and the fluctuations of the number of pairs approxi- 
mately 4 conforms the Poissonian law: 



w P {N C c) 



-(N cS ) (N C c) 



N r . 



N cE \ 



(61) 



Here w(N cS ) is the probability to observe N c£ charm 
quark- antiquark pairs in an event, provided that the av- 
erage number of cc pairs in this type of events is (iV cg ). 

The average number of particles with hidden, single, 
double, and triple charm is found by the convolution 



4 This approximation obviously breaks down at very large N C c 
when the energy of the produced cc pairs becomes comparable 
with the total energy of the system. But the probability of such 
events is clearly negligible. Only a tiny fraction of the total 
energy of the system is accumulated into the charm particles. 



of the results of the last section ({54"|) . (p l >?|) ~ ({I)9"jl with the 
probability (|61|) : 



(N k ) P = ™p(Ncc)(N k ) fl . 

N cS =l 

k = H, 1,1, 2, 2, 3, 3 



(62) 



(The subscript "P" stands for "Poisson"). Here we again 
see that the most of the charm is accommodated into the 
single charm hadrons: 



<A^ 0) - <A^ 0) 



(N c5 ). 



(63) 



The rest (a tiny fraction) is distributed over hidden, dou- 
ble and triple charm particles whose number can be easily 
found: 



(N c , 



{N c , 



' Nf 



(64) 
(65) 
(66) 



Surprisingly, the formulas for the number of double and 
triple charm particles l|65|l . (|66|l appeared to be exactly 
the same as in the grand canonical approach (|22[1 . (|23|> . 

However, this is not the case for the hidden charm. If 
the average number of cc pairs is small, (N C c) ^ 1, the 
average number of produced hidden charm hadrons 1|64|) 
is essentially larger than it would be naively expected 
from the grand canonical formula (|21|l . The two formulas 
(|21|l and (|64|l give similar results only when the number 
cc pairs is large, {N c£ ) 3> 1. The reason for this similarity 
is that the Poissonian distribution becomes narrow at 
{N c5 ) » 1: ((N c£ - {N c5 )) 2 ) = {N cS } « {N cz ) 2 . 

It is easy to see that any narrow probability distri- 
bution of the number of charm quarks and antiquarks 
would give the same result as the grand canonical ap- 
proach, provided that (N cS ) ^> 1. 



V. CHARM COALESCENCE IN A SUBSYSTEM 

The formulas of the previous section allow to calculate 
the number of charm particles in the entire system. It 
often happens, however, that only a part of the phase 
space is observed in the experiment: a limited rapidity 
interval, for instance. Formulas for charm production in a 
subsystem are necessary in this case. The charm and non- 
charm hadrons are distributed inhomogeneously in the 
(phase) space. Therefore, the number of observed charm 
particles in a subsystem cannot be calculated merely as a 
fraction of their total number in the system, proportional 
to the ratio of the volume of the subsystem to the total 
volume of the system. Moreover, the total volume of the 
system may be even unknown. The formulas derived in 
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this section are based on the thermodynamic parameters 
of the subsystem without any reference to those of the 
entire system. 

Let £ < 1 is the probability to find a charm quark in the 
subsystem, provided that exactly one cc pair is present 
in the entire system. Then, if the number of cc pairs in 
the entire system is iV cE , the probability to find N c charm 
quarks in the subsystem is given by the binomial law: 



w(N c \N cS ) 



N C \(N C£ -N C )\ 



£ N <>(1 (67) 



It is assumed that the distributions of quarks and anti- 
quarks are uncorrelated, and the probability distribution 
of the number of antiquarks N s is given by the same bi- 
nomial law. 

Then, if the total number of cc pairs in the system is 
fixed, the average number of the hidden charm particles 
produced in the system is given by a convolution of two 
probability distributions l|67[l with the right-hand side of 
JS3: 



(N H )fl = £ w(N c \N c5 ) w(N E \N cE )N c N s ^l- 



N c =0 N c =0 

Nh 



2 / at \2 



NiNj 



(68) 



In a more realistic situation, if the total number of 
cc pairs in the entire system fluctuates according to the 
Poissonian law H61[). the average number of hidden charm 
in the subsystem is given by 

(N H )T b = e(N cE )((N c5 ) + 1)-J^. (69) 
Similarly for the double and triple charm: 



<* 9 >a = m c -> 2 §, 

h 
Nf 



(70) 
(71) 



Note, that only the average number of the charm quark- 
antiquark pairs refers to the entire system in the above 
equations. The thermal quantities Nh, Ni and iVj are 
related to the subsystem under consideration. Therefore, 
one can calculate the number of charm particles produced 
in the subsystem, even if the thermodynamic parameters 
of the entire system are not known. Moreover, to calcu- 
late the number of double and triple charm, one does not 
actually need the total number of cc pairs in the entire 
system. It suffices to know this number for the subsys- 
tem £(-/V ce ). Only the calculation of the hidden charm 
requires the total number of charm pairs in the entire 
system, if this number is comparable to 1 or smaller. 



VI. THE CANONICAL APPROACH AND WHY 
IT IS NOT APPROPRIATE FOR THE CHARM 
COALESCENCE 

The canonical ensemble (CE) approach was initially 
proposed for a thermal treatment of the strangeness pro- 
ductions, when the average number of the strange parti- 
cles in the system is small (< 1) and the exact conser- 
vation of the net strangeness becomes important |l7j . It 
was also applied to the baryonic charge [lg ■ 

Although, as it will be explained later, this approach is 
not appropriate for the charm coalescence, I consider it in 
details because of two reasons. First, it has been widely 
used for the description of the charmonium production 
00 Second, it would be instructive to see, what 

changes, if the event-by-event fluctuations of the number 
of charm pairs is different from the Poissonian one. 

From the formal point of view, there is no problem to 
apply the canonical approach to the charm coalescence. 
Let us choose the probability distribution of the event- 
by-event fluctuations of the number of cc pairs as 



wce(N c5 ) 



j*Z NcE (V,T,{\}) 



where 



Z NcE (V, T, {A}) = Z NeN , (V, T, {A}) 



N c —Nc=N c 



(72) 



(73) 



In the zero approximation (|53[1 . the sum in the denomi- 
nator can be expressed via the modified Bessel function 



E ' 

iVcc=0 



2N cE 7 



(N c5 \y 



E 

JV cE =0 



I 2 7c VJViJV: 



(74) 



Again, to find the number of charm particles with dif- 
ferent charm content we have to convolute the probability 
((72^1 with the expressions (|S"lJl and (f57|) - lf!)9"|) similarly to 
(|62|l . The zero approximation results are expressed via 
the modified Bessel functions k — 0, . . . , 3: 



(Ni) { c ] E « (N-x)^ « (N c5 ) 



(75) 



h ( 2 lc a 'JN 1 N l 



I ( 2 7 fVJVi^I 



(76) 



? at J 2 ( ^x/NtNi 
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3 h ( 2 7 | 0) ViViiV i 



7 ( 2 7 i° ) V iV i iV i 



(78) 

The formula (|76[) for the hidden charm in CE has the 
same form as the corresponding formula in GCE ap- 
proach (see the middle part of J2U), but the result is 

different, because the value of 7 C °' > is not the same. In- 
deed, this value has to be found at given N cE from the 
transcendental equation (|75|l . Its solution may be quite 
different from (|T%jl . 

The equation (|75|l can be solved analytically in two 
limiting cases: (N cS ) ^S> 1 and (-/V ce ) <C 1. In the first 
case, the ratios of the Bessel functions tends to 1: 



h(x) 
Io{x) 



1 at x > I, k = 1, 2, 



(79) 



and the formulas are reduced to those of the grand canon- 
ical approach. 

In the second case, (N c£ ) -C 1, the Bessel functions can 
be replaced by the leading terms of their Taylor expan- 
sions: 



2 k k\ 



at x -C 1, k = 0,2, . 



The equation (|75|l then becomes 

(iv cc -)~ ( 7 (°)) 2 iv 1 iv I . 



The factor 7c can be easily found: 

v(°) 



7c 



' {N c5 ) 



(80) 



(81) 



(82) 



Then the number of hidden, double and triple charm par- 
ticles can be calculated: 



(Nh, 



(N 2 ) 



(N 3 ) 



(0) ~ 

CE — 




(83) 


(°) „ 
CE — 




(84) 


(0) „ 
CE — 


6 X cc/ Nl 


(85) 



The result (|83|l for the hidden charm coincides with 
that for the Poissonian fluctuation case (j64(l in the limit 
(N C c) <C 1. The reason for this coincidence is the follow- 
ing property of the both distributions: 



w{n) -C w(l) at (N ce ) < 1, n = 2, 3, 



(86) 



Therefore, the hidden charm hadrons are mostly pro- 
duced in the systems containing a single cc pair. The 



probabilities to observe exactly one cc pair becomes ap- 
proximately equal to each other for such distributions, if 
the average number of the pairs is the same and is small: 



w P (l) ~ wcbO) - (Ncc) at (N c5 ) 



e) < 1. 



(87) 



In contrast, the double and triple charm particles can- 
not be produced in the system containing only one cc 
pair. Two and three pairs at least are needed to pro- 
duce, respectively, a double and a triple charm particle. 
The probability to have two or three pairs in the system 
are different in the Poissonian and the canonical cases. 
From this reason, the results are essentially different. 

Although the Poissonian and canonical probability 
laws give the same result for the hidden charm in two 
limiting cases: (N cS ) -C 1 and (N cS ) 3> 1, they are up to 
about 10% different in the intermediate region (N c5 ) ~ 1 
□ 

Which of two approaches, the Poissonian or the Canon- 
ical one, is appropriate for the description of the charm 
coalescence in a system containing a small number of cc 
pairs? In fact, the key assumption (|72|l of the canonical 
approach relates the fluctuations of the number of charm 
pairs to the thermodynamic parameters of the system at 
chemical freeze-out. This is in an obvious contradiction 
with the basic postulates of the statistical coalescence 
model. Indeed, we have postulated that charm quarks 
are produced exclusively at the initial stage of the reac- 
tion in hard parton collisions. Creation and annihilation 
at later stages are neglected. Therefore, the fluctuations 
of the number of the cc pairs cannot have any relation to 
the properties of the system at the thermal stage. From 
these reasons, the canonical approach is not appropriate 
for the treatment of charm coalescence. 

The case of strangeness hadronization is essentially dif- 
ferent: strange quarks and antiquarks can be produced 
at the thermal stage, as far as the temperature is larger 
or comparable to the strange quark mass. Therefore, in 
the case of full strangeness thermalization (if not only the 
momenta, but the number of strange quark pairs is ther- 
mal), the canonical approach is the most appropriate. 
Still, it should be used with care if the full strangeness 
thermalization is not reached: j s ^ 1. 



VII. SUMMARY AND OUTLOOK 

The production of particles with single, double, triple 
and hidden charm in the framework of the statistical co- 
alescence model has been considered. The grand canoni- 
cal approach (Section^ is appropriate for systems con- 
taining a large number of charm quark- antiquark pairs: 
(N cS ) ^> 1- The solution can be found analytically l|20|) - 
(|23|l . if the number of hidden, double and triple charm 
particles is small comparing to the total charm. This 
is the case at presently available collision energies. At 
higher energies, the result can be found numerically from 
the coupled equations i|ll|) and (|12|l . 
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The grand canonical approach cannot be applied to 
the system containing a small number of charm quark- 
antiquark pairs: (N c s) < 1. It has been shown that the 
canonical approach is also inappropriate in this case. It 
is in variance with the basic postulates of the statistical 
coalescence model. As far as cc pairs are created in mu- 
tually independent nucleon-nucleon collisions, the fluc- 
tuations of N cS follows the Poissonian law. In this case, 
the result coincides with the grand canonical one for the 
double and triple charm |JBSJ|,(JBSJ, but differs essentially 
for the number of the hidden charm (|64ll . 

The charm coalescence in a part of a large system has 
been also studied. It is sufficient to know the thermal pa- 
rameters of the subsystem under consideration to find the 
number of particles with different charm content. The 
thermal properties of the entire system are not needed. 
Still one has to know the total number of the charm pairs 
in the entire system to calculate the number of hidden 
charm particles in the subsystem. For double and triple 
charm, the information on the number of cc pairs in the 
subsystem suffices. 

The obtained derived allow to obtain predictions for 



double, triple and hidden charm production in heavy ion 
collisions at all collision energies. At very high energies, 
like those of the Large Hadron Collider (LHC) in CERN 
(Switzerland), the expected number of produced charm 
quark-antiquarks is large (of order of hundreds or even 
more). Under these conditions, a sizable number of dou- 
ble and triple charm is expected. On the other hand, it 
would be interesting to see, whether the statistical coa- 
lescence model works at low energies. The possibility to 
study this at the accelerator facility with a very high lu- 
minosity which is planned to be build in GSI (Germany) 
is worth to be checked. A thermal or nonthermal behav- 
ior of heavy quarks can tell us much about the properties 
of the medium created during the heavy ion collision [2(| . 
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